Prom Neel long-range order to spin-liquids in the multiple-spin exchange model 
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The phase diagram of the multiple-spin exchange model on 
the triangular lattice is studied using exact diagonalizations. 
The two-spin (J2) and four-spin (J4) exchanges have been 
taken into account for 12, 16, 19, 21, 24, and 27 site samples 
in the parameter region J4 = — 0.25 (for a fixed J 2 = 1). It 
is found that the three-sublattice Neel ordered state built up 
by the pure two-spin exchange can be destroyed by the four- 
spin exchange, forming a spin-liquid state. The different data 
suggest that the phase diagram in this range of parameters 
exhibits two phases. The pure J2 phase is a three-sublattice 
Neel ordered phase, a small J4 drives it into a spin-liquid 
state with a spin gap filled of a large number of singlets. This 
spin-liquid phase is not of the same generic kind as the phase 
studied by Misguich et al. [Phys. Rev. B 60, 1064 (1999)]. 
It is observed on the finite-size samples that the spin liquid 
phase, as the Neel ordered phase, exhibits a magnetization 
plateau at m — 1/3, and for J4 > 0.15 a second plateau at 
m = 1/2. These two plateaus are associated respectively to 
the semi-classical orderings uud and uuud. 



I. INTRODUCTION 

The two-dimensional triangular lattice antiferromag- 
net (2D-TLA) was firstly proposed to be a candidate 
for the disordered (or spin-liquid) ground-state of the 
spin-i Heisenberg model (Anderson et al. in 70 's ||[|)- 
Different approaches failed to support this conjecture, 
but favor a ground-state with Neel long range order 
(LRO) • Nevertheless, the lattice frustration on 
the 2D-TLA attracts great interest of theorists, provid- 
ing a challenge for exotic antiferromagnets. Recently, 
the multiple-spin exchange model has been extensively 
studied as an alternative to the Heisenberg model, show- 
ing a rich structure of ground-states ]l|,|lO|,p]| . In this 
model, the ground-state can be ferromagnetic (FM), anti- 
ferromagnetic (AFM) with Neel LRO, or a spin- liquid 
(SL). A prospective phase diagram has been given by 
Misguich et al., who considered two-, four-, and five-spin 
exchange interactions on the 2D-TLA |Q. They found 
that a large enough four-spin exchange interaction drives 
th FM phase into a SL phase. They did not study how 
the AFM Neel LRO is destroyed by the four-spin ex- 
change interaction, and how the transition between Neel 
LRO and the short range RVB phase takes place. This 
question is the main object of this paper. 



Unhappily there is no exact method allowing the study 
of the zero temperature phases of such frustrated sys- 
tems. For a finite system, however, one can always, in 
principle, represent the eigenstates in the complete ba- 
sis of spin configurations. This allows one to have a real 
touch on the exact ground-states of small size systems 
through numerical computations. The huge number of 
spin configurations (2 N ) becomes a great obstacle on the 
way of numerical simulations. On the most recent com- 
puters, the largest sample that may be handled in exact 
diagonalizations has 6x6 sites. On the triangular lat- 
tice, the Quantum Monte Carlo method is plagued by 
the well-known sign problem, but a new technique called 
Stochastic Reconfiguration allows handling samples up 
to 12 x 12 sites All these calculations point to Neel 
LRO, with a sublattice magnetization of the order of 40% 
of the saturated value ||. Series expansions give a re- 
duced (20%) but non-zero sublattice magnetization Q . 

On the other hand, in the case of short range correla- 
tions, the situation is more straightforward, as soon as 
the available sizes are of the order of, or larger than the 
correlation length. This is fortunately the case in the SL 
phase found in the J2 — </ 4 model ( J2 < 0, J4 > 0) by 
Misguich et al. gQ. 

In this work, we use exact diagonalizations to obtain 
the exact eigenenergies versus wave vectors and total spin 
for 12, 16, 19, 21, 24 and 27 site samples of the J 2 - J 4 
model ( J2 = 1, J4 > 0) on the 2D-TLA. In the classical 
limit, the AFM ground-state of the 2D-TLA can be de- 
scribed as a three-sublattice structure, with spins of dif- 
ferent sublattices making angles of 27r/3. Periodic bound- 
ary conditions are compatible with the three-sublattice 
structure for samples with 12, 21, 24 and 27 sites, but not 
for the 16 and 19 site samples. Therefore, we use twisted 
boundary conditions for the 16 and 19 site samples ||] 
and periodic boundary conditions for the 12, 21, 24 and 
27 site ones. 



II. THE MODEL: THE MULTIPLE-SPIN 
EXCHANGE HAMILTONIAN 

The Hamiltonian of the multiple-spin exchange model 
is given by 

H = E(- 1 )' lJ «( P « + P n ') . J n > ,n > 2 (1) 
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where J„ are the n-spin exchange tunneling probabili- 
ties (exchange coefficients), P n and P~ x are the n-spin 
exchange operators and their inverse operators, respec- 
tively. The alternative sign in the summation over n in 
Eq. [I] comes from the permutation of fermions. In gen- 
eral, the exchange coefficients decrease with increasing n. 
The two-spin exchange term gives exactly the Heisenberg 
Hamiltonian up to a constant, since one has 



Pi = 2S, ■ S,- 



(2) 



where Si and s^- are spins localized at site i and j, re- 
spectively. The three-spin exchange operator is exactly 
equivalent to a sum of two-spin exchange operators |Q. 
Thus, the three-spin exchange term of equation ([!]) can 
be absorbed into the two-spin exchange term, as long as 
J2 is replaced by the effective two-spin exchange coeffi- 
cient J£ ff — Ji~ 2J 3 . Therefore, except for the two-spin 
exchange, the next most important term is the four-spin 
exchange. A pure positive two-spin exchange (i.e., the 
Heisenberg Hamiltonian) on the 2D-TLA gives an AFM 
phase with Neel LRO, and, as shown in Ref. [Q, a pure 
four-spin exchange gives a SL phase. In this paper we 
use the specific properties of the spectra of these different 
kinds of phases to study the transition from one phase to 
the other, when the relative weight of the four-Spin Ex- 
change J4 increases relatively to the antiferromagnetic 
two-spin coupling. In the following, all the energies are 
measured in units of J| ff = 1. 

Little is known on this region of the phase diagram. 
Previous works are based on a classical approximation |Q , 
semi-classical spin-wave calculations [|| or mean-field 
Schwinger-boson results 0. The classical result pre- 
dicts a transition from the 3-sublattice Neel state to 
a 4-sublattice tetrahedral state at J4 = 0.24. Both 
quantum approaches indicate that the four-spin exchange 
strongly enhance fluctuations in the 3-sublattice Neel 
phase. Kubo et. al [0 found that the sublattice mag- 
netization vanishes for J4 > 0.17. In the Schwinger- 
boson approach j3|, the Neel state is destroyed when 
J4 > 0.25. These two techniques have a general tendency 
to underestimate the effects of quantum fluctuations on 
ordered phases. The exact diagonalization analysis pre- 
sented here indeed shows that the Neel long-ranged order 
disappears for a smaller value of J4 (the critical value is 
estimated to be in the interval Jf 7 ~ 0.07 • • ■ 0.1). 



III. CRITERION TO DISCRIMINATE BETWEEN 
NEEL LRO AND SL PHASE: THE SPIN GAP ? 

A. Finite-size energy spectrum of the Neel LRO 
phase 

In the classical limit, a A-site 2D-TLA sample with 
Neel LRO is characterized by a three-sublattice structure 



with spin A/6 on each sublattice. Coupling of these three 
A/6-spins, gives total spin S with min{2S +1, N/2 — S + 
1} degeneracy [||. In an isotropic antiferromagnet (as the 
collinear AFM which has equal spin susceptibilities and 
spin wave velocities) the finite-size total energy depends 
on the total spin 5 (to first order in 1/A) as: 



E S = E Q 



2N X 



5(5+1), 



(3) 



where Eq = Aeo is the energy of the ground-state in 
the thermodynamic limit and x is the isotropic magnetic 
susceptibility of the sample. In the anisotropic case this 
equation should be rewritten: 



Es — En 



i^S(S+l) + ^-(--—)S!, (4) 

2Ax_l 2A xil X-L 



where S3 is the component of the total spin 5 on the in- 
ternal symmetry axis of the spin system, and x\\ an d Xj_ 
are the magnetic susceptibilities on the internal symme- 
try axis and on the perpendicular plane, respectively. In 
the broken symmetry picture the symmetry axis is per- 
pendicular to the plane of the spins and \\\ (respectively 
X_i_) measures the spin fluctuations orthogonal to (respec- 
tively in) the spin plane. Eq. || (respectively Eq. p, is 
the dynamical equation of a rigid rotator (respectively of 
a quantum top). 

Eqs. P show that the slopes of the total energy 
versus 5(5+1) and 5| approach zero as 1/A does when 
A — > 00. 53 is an internal quantum number dynami- 
cally generated, which is not under control in a finite-size 
study. But the total spin 5 is a good quantum number 
and the A -1 scaling of the 5(5 + 1) dependence of the 
total energy versus sample size is interesting because it is 
more rapid than the scaling law of the order parameter 
(which goes as A -1 / 2 ) IMQ. 



B. Finite-size scaling in the SL phase 

In a SL phase, contrarily to the Neel LRO phase, the 
spin gap (i.e. the difference in total energy between 
ground-states in the 5 = 1 sector and in the 5 = 
sector) does not collapse to zero in the thermodynamic 
limit. The finite-size scaling law in this second situation 
is not known exactly, insofar as the "massive" phase is 
not characterized precisely. Heuristically, we expect the 
finite-size spin gap to decrease exponentially to a finite 
value A (00) with the characteristic length £ of the spin- 
spin correlations. For samples of linear size L smaller 
than the correlation length and in the cross-over regime, 
there are not enough quantum fluctuations to destroy 
the sublattice magnetization and the system probably 
behaves as it were classical (i.e. with a spin gap decreas- 
ing as A -1 ). The following heuristic law might be used 
to interpolate between the two behaviors: 
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A(L) = A(oo) + ^ x exp{-L/i) 



(5) 



IV. SYMMETRIES OF THE LOW LYING 
LEVELS IN A NEEL ORDERED PHASE 



C. Quantum critical regime 

The use of this heuristic law (Eq. |5|) encounters a severe 
difficulty as soon as the disordered system approaches a 
quantum critical point: in such a situation the correlation 
length £ diverges, the gap closes to zero and on a finite- 
size sample it is impossible to discriminate between such 
a situation and isotropic Neel LRO (Eq. ||) . 

D. Numerical results 

In view of this difficulty we have done a pedestrian 
finite-size scaling of the spin gap by using the simplest 
linear 1/N behavior which probably gives a lower bound 
of the gap in a SL outside critical points. The physical 
reason is that we expect the finite-size corrections of the 
gap value to be smaller in a system with a finite corre- 
lation length (no long-range order) than in a LRO Neel 
phase where it vanishes as 1/N. [] The results extrap- 
olated to N — > co are shown in Fig. [I]. Strictly speak- 
ing the spin gap never extrapolates to zero except for a 
pure Ji where it is equal to zero within its error bar (At 
J4 = 0, data of Bernu et al. |q| for TV = 36 are added to 
the present results, see Fig. fiT). Nevertheless these data 
already show three distinct ranges for the parameter J4: 
for very small J4 (below 0.075) Neel LRO is plausible but 
should be confirmed by another approach. For J4 larger 
than 0.1 a gap certainly opens rapidly with increasing J4 
and then decrease for J4 > 0.175. The spin gap criterion 
cannot give more insight on the phase diagram. We will 
now move to the analysis of the symmetries of the low 
lying levels of the spectra to characterize more precisely 
these three phases. 



1 This assumption would be invalid at a critical point where 
the uniform susceptibility vanishes: In such a case the gap 
might close as 1/ y/N. The analysis of section |^ shows that 
the behavior of the system changes rather abruptly from a 
Neel like spectrum to a spectrum with a very large number 
of low lying singlets below the first S — 1 state (and poten- 
tially a T = residual entropy in the singlet sector). This 
happens before any decrease to zero of the spin velocity or of 
the homogeneous spin susceptibility. So, if the true thermo- 
dynamical system has indeed a critical point between the two 
phases, the sizes we are looking at are too small to scrutinize 
the critical regime. 



A. Theoretical background 

Firstly we show the low energy spectrum of the pure 
Heisenberg model on the 21 site sample (Fig. |J). In order 
to emphasize the low energy structure we have displayed 
the low energy spectrum minus a rigid rotator energy 
aS(S+l). 

Let us first concentrate on the lowest part of the energy 
spectrum in each S sector (solid and open triangular sym- 
bols in the figure). This family of levels forms on a finite- 
size lattice the quantum counterpart of the semi-classical 
Neel state. These specific states are in the trivial repre- 
sentation of the invariance group of the three-sublattice 
Neel ordered solution j8|. To be definite: 

• The Neel ground-state breaks the I-step translation 
but is invariant in a 3-step translation: as a con- 
sequence the only wave vectors appearing in this 
family of QDJS (for Quasi-degenerate-joint-states 
defined by Bernu et al. ||) are respectively the cen- 
ter k = (0, 0) and corners ±ko of the Brillouin zone. 

• These QDJS belong specifically to the trivial rep- 
resentation of Cz v as the Neel state itself (i.e. they 
are invariant in a 2tt/S rotation, and in a reflection 
symmetry). 

• As the 7T rotation symmetry of the lattice is broken 
in this particular ground-state, these QDJS appear 
cither in the odd or even representation of the 2- 
fold rotation group (see Ref. M for more details). 

• The numbers and characteristics (quantum num- 
bers) of the QDJS in each S sector are precisely 
fixed by theory M: for the 21 sites spectrum dis- 
played in Fig. g(as for all sizes that have been 
studied up to now) the numbers of low lying levels 
and their quantum numbers correspond exactly to 
the above-mentioned theoretical predictions. 

B. Numerical results 

The dynamical law given by Eq. |] is still imperfectly 
obeyed for the 21 site sample: in particular the genera- 
tion of the internal symmetry is still imperfect but nev- 
ertheless the spectrum of a quantum top could already 
be anticipated. 

Above these levels with specific properties, there ap- 
pear eigenstates with wave vectors belonging to the inside 
of the Brillouin zone (simple dashes in Fig. |J). A group 
of such eigenstates with different total spin represents 
a magnon excitation of the Neel ground-state. As the 
antiferromagnetic magnons have a linear dispersion law, 
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the softest magnon energy scales as the smallest wave 
vector accommodated in the Brillouin zone of the finite- 
size sample. Thus, for sizes large enough, these levels 
collapse to the ground-state as 1/ y/N, more slowly than 
the QDJS which collapse as 1/N to the thermodynamic 
Neel ground-state energy. That is the reason of the ap- 
pearance in Fig. |^ of a quasi gap between the QDJS and 
the magnon excitations. 

This hierarchy of low lying levels is a very strong con- 
straint on the finite-size samples spectra. It is perfect for 
sizes up to 27 and for J4 smaller or equal to 0.075 and 
totally absent for J4 larger or equal to 0.1. 

• This result, associated to the spin gap behavior, 
consistently proves that for J4 larger or equal to 
0.1 the system is in a spin-liquid state with rather 
short range spin-spin correlations. 

• For J4 < 0.075 the structure of the low lying eigen- 
levels of the spectra are compatible with Neel LRO. 
BUT as discussed above, it is indeed impossible to 
precisely point a quantum critical transition within 
this approach. In view of the spectra, we might 
speculate that the transition is second order and 
that it takes place between 0.07 and 0.1. 

• From J4 = to J a ~ 0.1 we see a softening of the 
spin-wave velocity consistent with the gradual de- 
crease of the Neel LRO (Fig. ||). However, Fig. || 
shows that a large number of singlet states are al- 
ready present at low energy when J a = 0.1. There- 
fore, the system is certainly no longer in the or- 
dered phase at J4 ~ 0.1. We can conclude from 
Fig. [| that the spin wave velocity does not vanish 
at the critical point (even if the precise location 
of the transition cannot be determined). This has 
been previously suggested by various analytical ap- 
proaches |15|-|l7|| . The sizes studied are neverthe- 
less too small to check Azaria's prediction |l5|] of 
an 0(4) symmetry of the effective field theory at 
the critical point. 

V. THE LOW ENERGY EXCITATIONS OF THE 
SL PHASE 

Contrarily to our expectations, the SL phase which 
appears immediately after the disappearance of the Neel 
ordered phase is not the phase studied by Misguich et 
al. 0. It is indeed a SL phase, with a gap and short 



range spin-spin correlations. But, as might be seen in 
Fig. |], this phase exhibits a very large number of singlets 
in the magnetic gap and seems in this respect similar 
to the spin-liquid phase of the Heisenberg model on the 
kagome lattice |l8|,|l9|]. However, we are not aware of any 
exponential degeneracy (i.e. ~ exp(iV)) in the classical 
MSE model [J as is the case for the classical kagome 
antifcrromagnet . 

We suspect that a much larger four-spin exchange pa- 
rameter will be needed to recover the SL phase stud- 
ied by Misguich et al. These data point to the exis- 
tence of this new phase in a finite range of parameters 
0.075 < J4 < 0.25. However, one cannot disregard the 
hypothesis that these properties are in fact those of a 
critical point, with a critical region enlarged by finite- 
size effects. More work with different methods is needed 
to clarify this point. 

VI. MAGNETIZATION PLATEAUS 

In an external magnetic field B along the z axis, the 
total energy of the state with component S z of the total 
spin is given by: 

E B = E S - S Z B (6) 

The magnetization is determined by the minimum of Eb 
respective to S z , which requires dEs/dSz = 0. There- 
fore, in the isotropic case with Neel LRO, one has 

to = 2x5 - 1/N (7) 

where m = 2S/N, is the polarization relative to the sat- 
urated magnetization N/2 and x is indeed the magnetic 
susceptibility. 

We noticed that when the four-spin exchange interac- 
tion increases, deviation from Eq. |^ occurs about S z — 
N/6 and N/4. This is in agreement with the earlier 
mean-field calculation of Kubo and Momoi || who pre- 
dicted magnetization plateaus at m = 1/3 and m = 1/2 
in the J2 — J a model. We present the magnetization 
curve of the 24-site sample in Fig. ||. A small plateau 
exists at 1/3 magnetization for J4 = 0, and its width 
first increases slowly as J a increases and then decreases 
from around J4 = 0.125. This 1/3 plateau has also been 
found in previous studies (see [ElJ and references therein) 
in the pure three-sublattice Neel ordered system. The 
m = 1/2 magnetization appears at about J4 = 0.1 and 



2 An illustration of the striking qualitative differences be- 
tween spectra in the Neel region and in the spin-liquid phase 
has been previously published: two N = 27 spectra at J4 = 
and J a — 0.1 are displayed in Fig. 2 and Fig. 3 of Ref. |j. 



3 Momoi et al. |2(J found a ground-state degeneracy in the 
classical MSE model but this degeneracy only grows as the ex- 
ponential of linear size of the system (~ exp s/N). Moreover, 
this classical degeneracy was only found in a region where 
J-2 < is ferromagnetic (| < J 4 /| J 2 | < |). 
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its width increases with J4. Finally, there still exists a 
plateau at about m = 1/2 in samples with odd numbers 
of sites, but it distributes over the two closest positions 
to the 1/2 magnetization. The m = 1/3 and m = 1/2 
plateaus correspond to the classical uud and uuud order- 
ing of spins (see Momoi et al. |2(J for the four-sublattice 
uuud state). These phases are more "classical" than the 
zero-field phase, but yet show a decrease of the sublattice 
magnetization from the classical saturation values. 



VII. CONCLUSION 

In conclusion, we studied the transition between the 
Neel ordered state and a Spin Liquid state of the 
multiple-spin exchange model by means of the exact di- 
agonalization method. The pure three-sublattice Neel 
ordered phase is gradually destroyed by quantum fluc- 
tuations when increasing the 4-spin exchange coupling. 
The spin wave velocity decreases but apparently remains 
finite at the transition. The quantum disordered phase 
tuned by the 4-spin exchange coupling is different from 
the pure J4 phase studied by Misguich et al. It exhibits 
low energy singlet excitations, reminding of kagome SL. 
This result opens many interesting questions that cannot 
be answered in the present framework: is this phase a new 
generic SL phase or a finite-size manifestation of a quan- 
tum critical regime ? Are these singlet excitations the 
"resonon" modes invoked by Rokhsar and Kivelson |^2j ? 
In agreement with previous studies pi , we find two mag- 
netization plateaus at 1/3 and 1/2 of the full magnetiza- 
tion. These plateaus are associated to the semi-classical 
uud and uuud ordering structures. A finite-size scaling 
on much larger sizes is needed to draw a definite conclu- 
sion on this magnetic phase diagram. 
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FIG. 1. Extrapolation N — > 00 of the spin gaps of 12, 21, 
24 and 27 site samples for different J4. The inset shows the 
extrapolation at J4 = 0. The spin gap for N = 36 is evaluated 
from data of Ref. [B. Since not all symmetry sectors were 
investigated for S = 1 and N = 36, the energy of the lowest 
triplet is not known exactly. Therefore we estimated the spin 
gap by A ~ i [E(S = 3) - E(S = 0)], as in Ref. §. Lines are 
guides for the eye. 
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FIG. 2. Spectrum of the Heisenberg model ( J a = 0)for the 
21-site sample. v : levels with k = (0,0) and symmetries 
^2ir/3 = 1 an d Rw = 1 {Re is the phase factor obtained 
in a ^-rotation about the origin); black solid v : levels with 
k = (0,0) and symmetries R^-n/s = 1, Rn = — 1; >: levels 
with ko (the corners of the Brillouin zone) and symmetries 
^2ir/3 = 1; — : levels with wave vectors inside the Brillouin 
zone. 



FIG. 3. Spin wave velocities of 21, 24 and 27 -sites samples. 
Lines are guides for the eye. These quantities are computed 
as the ratio of the first AS = 1 excitation energy divided by 
the momentum of the corresponding excitation (according to 
v = A-E/Ak). For Ji > 0.1, these numbers do not correspond 
to well defined physical excitations because of the strong per- 
turbations of the spectra. Notice at J4 = 0.1, Fig. ^ exhibits 
a huge number of low-energy singlet states below the spin gap 
and the system is therefore already in a SL phase. 
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FIG. 5. Magnetization versus external magnetic field for 
different J 4 (N = 24). The dotted lines label the 1/3 and 1/2 
magnetization. The top panel shows the widths of the 1/3 
(V) and 1/2 (□) plateaus versus J a 
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